Multi-objective exponential distribution optimizer (MOEDO): a novel math-inspired multi-objective algorithm for global optimization and real-world engineering design problems

The exponential distribution optimizer (EDO) represents a heuristic approach, capitalizing on exponential distribution theory to identify global solutions for complex optimization challenges. This study extends the EDO's applicability by introducing its multi-objective version, the multi-objective EDO (MOEDO), enhanced with elite non-dominated sorting and crowding distance mechanisms. An information feedback mechanism (IFM) is integrated into MOEDO, aiming to balance exploration and exploitation, thus improving convergence and mitigating the stagnation in local optima, a notable limitation in traditional approaches. Our research demonstrates MOEDO's superiority over renowned algorithms such as MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO. This is evident in 72.58% of test scenarios, utilizing performance metrics like GD, IGD, HV, SP, SD and RT across benchmark test collections (DTLZ, ZDT and various constraint problems) and five real-world engineering design challenges. The Wilcoxon Rank Sum Test (WRST) further confirms MOEDO as a competitive multi-objective optimization algorithm, particularly in scenarios where existing methods struggle with balancing diversity and convergence efficiency. MOEDO's robust performance, even in complex real-world applications, underscores its potential as an innovative solution in the optimization domain. The MOEDO source code is available at: https://github.com/kanak02/MOEDO.

(OAs) and introduce more potent OAs 3 .Thanks to technological progress, newer OAs that boast superior efficiency, precision and speed in addressing diverse optimization tasks are gaining traction [4][5][6] .Moreover, specific challenges like local optima and the irregularities and non-convexities of exploration domains have played a pivotal role in this evolution.
Such constraints on OAs have spurred scholars and industry experts to develop innovative metaheuristic algorithms to navigate diverse optimization hurdles 7,8 .Parallel to the growth of information technology, a plethora of optimization challenges have emerged across sectors like engineering, bioinformatics, operations research and geophysics 9,10 .Numerous optimization issues are categorized as NP-hard, implying their solutions are not achievable within polynomial time unless NP is equivalent to P. As a result, exact mathematical methods are typically reserved for problems of a smaller scale.
Researchers have explored alternative strategies (approximation techniques) to identify feasible solutions within a reasonable timeframe rather than abandoning the effort.These techniques can be broadly categorized into heuristics and metaheuristics.The primary difference between the two is that heuristics are closely tied to the specific nature of a problem, making them effective for particular challenges but less so for others.In contrast, metaheuristics present a more universal algorithmic structure or a black-box optimization tool suitable for almost any optimization problem (OP).These advanced heuristics, designed to address a range of OPs, are termed metaheuristics (MHs).Numerous metaheuristic algorithms (MHAs) have been successfully employed in recent times to navigate complex challenges 11 .A key benefit of these algorithms in addressing intricate optimization tasks is their capability to identify commendable solutions, irrespective of the problem's scale or intricacy 12 .MH algorithms have been utilized across a spectrum of optimization challenges, encompassing both single and multi-objective, as well as continuous and discrete scenarios 13 .
In the real world, problems often present multiple objectives, some of which might conflict with others.As a result, multi-objective optimization problems (MOPs) align more closely with these multifaceted challenges than single-objective optimization does 14 .A common approach to solving MOPs is to treat each objective as an individual OP, addressing them one after the other based on their respective significance 15 .Alternatively, individual solvers can exchange insights across successive iterations.
Techniques for multi-objective optimization can be divided into three primary categories 16,17 : a priori, a posteriori and interactive approaches.
• A priori techniques: these strategies rely on preliminary knowledge about the problem and its objectives.They aim to pinpoint the Pareto optimal solutions even before the optimization process kicks off.A prevalent tactic in these techniques is transforming a MOP into a single-objective challenge at the outset.Notable examples of a priori techniques encompass weighting methods, goal-oriented programming and lexicographic strategies 18 .• A posteriori technique: these strategies draw from the results secured during the optimization phase.They strive to determine the Pareto optimal solutions post the completion of the optimization.Evolutionary algorithms like genetic algorithms, particle swarm optimization and simulated annealing fall under this category, as do Multi-objective Genetic Algorithms such as NSGA-II, SPEA2, MOEA/D and NSGA-II 19 .• Interactive techniques: these strategies necessitate human engagement throughout the optimization phase.
They grant the decision-maker the ability to engage with the optimization process and offer insights on the outcomes produced by the algorithm 20 .Examples of interactive techniques include interactive genetic algorithms, interactive evolutionary tactics and interactive particle swarm optimization.
While a priori techniques are generally straightforward to deploy and can yield prompt results, they might miss out on capturing the genuine Pareto optimal solutions.On the other hand, a posteriori technique, though potentially offering more precise outcomes, can be intricate to implement and might be resource-intensive.Interactive techniques strike a balance, enabling the decision-maker to influence the outcomes, potentially enhancing solution quality.The techniques are illustrated in Fig. 1.
Multi-objective evolutionary algorithms (MOEAs) have gained popularity as effective a posteriori technique for addressing MOPs.Their strength lies in their ability to approximate the Pareto set (PS) and/or Pareto front (PF) within a single execution 21 .A standard MOEA operates in two main stages for every generation.The first stage employs evolutionary mechanisms like crossover and mutation to generate new solutions in the decision realm.The second stage involves selecting superior solutions from both the newly created and pre-existing solutions for the subsequent generation 22 .Depending on the selection methodologies employed, most contemporary MOEAs can be categorized into four types: domination-based, indicator-based, decomposition-based and hybrid MOEAs 23,24 .
To determine the quality of selection, many MOEAs assess the function values of all the freshly produced solutions.Consequently, these algorithms often necessitate numerous function evaluations to achieve convergence 25 .Yet, the long-term outcomes of traditional optimization methods leave much to be desired 26 .Issues such as the sensitivity to the problem's initial estimation, reliance on the precision of the differential equation solver's solution and the risk of becoming ensnared in local optima conventional optimization techniques 27 .In addressing highly nonlinear challenges, there is a heightened risk of settling for a solution that is merely locally optimal 28 .
Metaheuristics for MOPs, commonly referred to as MOEA, can be categorized based on their solution selection methodologies into four primary groups: dominance-based, decomposition-based, indicator-based and hybrid selection mechanisms.
1. Dominance-based MOEAs.These prioritize achieving optimal convergence by assigning fitness scores to individuals according to the Pareto-dominance principle.Additionally, a specific strategy is provided to maintain diversity.A widely recognized approach in this category ranks individuals through rapid nondominated sorting, emphasizing reduced time complexity [29][30][31][32] .When multiple individuals possess the same non-domination rank, density metrics ensure diversity, leading to a more uniform distribution of solutions.However, as the number of objectives increases, the proportion of non-dominated solutions grows exponentially, making Pareto-dominated MOEAs less suitable for many-objective optimization.This can hinder the evolutionary process, diminishing selection pressure, diversity and convergence.2. Decomposition-based MOEAs.These segment a MOP into multiple scalar optimization sub-tasks using a decomposition strategy.Solutions for each subtask are optimized by conducting evolutionary operations among its various neighboring sub-problems.The most frequently employed decomposition-based MOEA is MOEA/D, introduced by Zhang and Li 30 .The MOEA/D structure can incorporate various traditional single-objective optimization and localized search methods [33][34][35] .3. Indicator-based MOEAs.These MOEAs utilize performance metrics related to solution quality as selection standards, guiding the search towards consistent enhancement of the overall population's anticipated attributes 36,37 .IBEA, a pioneering indicator-based MOEA, was developed by Zitzler and Kunzl 38 .It employed a binary indicator to evaluate the comparative benefits of two approximation solution sets, aligning with the decision-maker's preference.Given its compatibility with the Pareto principle, this indicator can be employed to determine fitness.This approach offers a foundational model for future indicator-based MOEAs.Some studies suggest that HV indicators might be a more computationally efficient alternative to other indicators while retaining similar theoretical properties.4. Hybrid MOEAs.Addressing intricate MOPs often necessitates leveraging the strengths of multiple algorithms, leading to the introduction of various hybrid MO algorithms [39][40][41][42] .The primary considerations for hybrid algorithms revolve around the selection of algorithms to merge and the methodology for their integration.Contemporary research predominantly centers on amalgamating different strategies to generate offspring populations.However, there is a need for further exploration in other domains.
Several algorithms have been crafted to adeptly navigate intricate OPs.One notable example is the exponential distribution optimizer (EDO) algorithm introduced by Mohamed Abdel-Basset et al. 64 .This innovative algorithm has been applied across diverse sectors, from engineering challenges to manufacturing processes and scientific modeling.The foundational principle of the EDO algorithm draws integrates two approaches geared towards exploitation and exploration tactics.During the exploitation phase, the algorithm employs three core www.nature.com/scientificreports/principles: the absence of memory property, a guiding solution and the exponential disparity observed among the exponential random variables to refine the existing solutions.
In recent times, the exponential distribution optimizer (EDO), a math inspired metaheuristic method rooted in the fractal theorem, was presented.Our research aims to introduce a novel multi-objective iteration of the EDO algorithm tailored for MOPs.The proposed version integrates the NDS, CD and IFM principle and a stochastic learning mechanism.To gauge the potency of this proposed method, we employ distinct benchmark test functions: ZDT 65 , DTLZ 66 , Constraint 67,68 (CONSTR, TNK, SRN, BNH, OSY and KITA) and real-world engineering design Brushless DC wheel motor 69 (RWMOP1), Helical spring 68 (RWMOP2), Two-bar truss 68 (RWMOP3), Welded beam 70 (RWMOP4), Disk brake 71 (RWMOP5).The objective of this assessment is to compare the efficacy of our proposed method against MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO, using metrics like generational distance (GD) 34 , inverse generational distance (IGD) 35 , hypervolume 36 , Spacing 37 , Spread 36 and run time (RT).The approximations of the Pareto-front produced by our method are evaluated using these metrics.
The key contributions of this research are succinctly outlined as follows: 1. Introduction of the multi-objective exponential distribution optimizer (MOEDO) Algorithm, incorporating non-dominated sorting (NDS) and crowding distance (CD) principles.2. Integration of the information feedback mechanism (IFM) to decompose multi-objective challenges into single-objective sub-tasks, enhancing the algorithm's efficiency.3. Utilization of the IFM approach to ensure a balanced dynamic between exploration and exploitation, fostering improved convergence and the capability to bypass local minima.4. Comprehensive evaluation of MOEDO's performance against established multi-objective methods using benchmark datasets like ZDT, DTLZ, CONSTR and real-world engineering design problems.5. Employment of metrics such as generational distance (GD), inverse generational distance (IGD), hypervolume (HV), spacing, spread and run time (RT) to assess results, demonstrating the superior capabilities of our proposed technique in various test scenarios.
The structure of our comprehensive research is as follows: "Methodology" elaborates on the exponential distribution optimizer."Results and discussion" provides an in-depth look at our proposed MOEDO method for multi-objective global optimization."Conclusion" presents our experimental results.

Methodology
This section delves into the EDO algorithm 64 fundamental components.We will first discuss the inspiration behind EDO, followed by a detailed explanation of its initialization.We will also examine the explorative and exploitative features of EDO before discussing its mathematical structure.

Exponential distribution model
The EDO approach takes its foundational principles from the theories of exponential distribution.This distribution, which is continuous in nature, has been instrumental in explaining several occurrences in the natural world.For instance, the time span between the present moment and the occurrence of an earthquake can be represented using this distribution.Similarly, the time-based likelihood of a car reaching a toll station aligns with the exponential distribution.Exponential random variables have frequently been utilized to shed light on past events, particularly focusing on the duration leading up to a specific incident.Here, we delve into the mathematical construct of the exponential distribution and elucidate its unique attributes.Imagine having an exponential random variable, represented by x , associated with a parameter, .This relationship can be expressed as x ∼ EXP( ) .The Probability Density Function (PDF) of this variable is indicative of the duration by: As time is a continuous factor, it must always hold a non-negative value, i.e., (x ≥ 0) .Significantly, the parameter, which is always positive, denotes the rate of occurrence in the exponential distribution.The exponential distribution cumulative distribution function (CDF) can be derived through a specific formula: A higher value of the rate of occurrence, , indicates a reduced probability of the concerned random variable.CDF function ascends, commencing from the foundational exponential rate and intensifying in proportion to the growth of the exponential random variable.For an exponentially distributed random variable, its mean (µ) and variance σ 2 can be articulated using specific mathematical expressions: These expressions reveal an inverse relationship between the parameter and both the mean and variance.To put it succinctly, as increases, both the mean and variance decline.Moreover, the standard deviation (σ ) mirrors the mean value and its derivation follows a particular computation (1) Vol www.nature.com/scientificreports/ The memoryless nature of exponential distribution One of the unique characteristics of certain statistical probability distributions is the 'memoryless' attribute.This implies that the chance of a forthcoming event transpiring is independent of past events.In simpler terms, previous occurrences do not influence future probabilities.The exponential distribution, which is a continuous type, encapsulates this memoryless attribute, especially when gauging the timespan before an event occurrence.When a random variable, denoted by x , adheres to the exponential distribution with this memoryless feature, it means for any positive whole numbers s and t belonging to the series {0, 1, 2, . . ., ∞} the following holds true: Launching with the initial group During the onset or the initiation stage, we cultivate a group termed ( X winners ), comprising N diversely valued, randomly formed solutions.To depict this search operation, we utilize an assortment of exponential distributions.Every potential solution is perceived as an embodiment of the exponential distribution.The respective positions of each solution are considered as exponential random variables conforming to this distribution.These are then structured as vectors, having a dimension d.
In this setup, X winners symbolizes the j element of the ith candidate within the exponential distribution vector.Following this, we can define the preliminary group, X winners .
To randomly generate each variable within the candidate exponential distribution in the solution space, we utilize a particular formula.
Here, lb and ub delineate the lower and upper limits of the given problem.The term 'rand' signifies a randomly derived number within the span [0,1].Upon concluding the initiation, we embark on the optimization phase.This leans on the exploratory and refinement skills of our proposed approach over iterative rounds.Subsequent sections will elucidate the two core techniques (exploration and refinement within EDO) for pinpointing the global pinnacle.

EDO exploitation strength
The exploitation component of EDO harnesses various facets of the exponential distribution model, including its memoryless trait, exponential rate, typical variance and average value.Additionally, a guiding solution steers the search towards the global peak.Initially in EDO, a set of random solutions is cultivated, mimicking an array of exponential distribution patterns.These solutions undergo evaluation via an objective function and are subsequently ranked in terms of performance.For maximization problems, they are ordered in descending efficiency, while for minimization challenges, they are arranged in ascending order.Areas surrounding a robust solution are fertile grounds for pinpointing the global apex.This is why several algorithms probe spaces around robust solutions, drawing weaker ones towards them.Therefore, the global apex quest focuses on the guiding solution.The guiding solution, labeled ( Xguide time ), is derived from the mean of the top three solutions in a sorted group.
In this context, time replaces the term iteration, alluding to the period until the subsequent event in the exponential distribution.The maximum number of iterations is symbolized by Max_time.This guiding solution offers invaluable insights about the global apex.Instead of exclusively leaning on the current best solution, the guiding one is prioritized.Though the leading solution has vital details about the global peak, an exclusive focus on it could lead to convergence around a local maximum.By introducing a current guiding solution, this challenge is mitigated.To emulate the memoryless characteristic of the exponential distribution, a matrix, termed memoryless, is formulated.This matrix houses the latest solutions, irrespective of their present efficiency scores.Initially, it mirrors the original population.Following this, solutions generated at the present time, regardless of their efficiency, are stored in the matrix, dismissing their past contributions.Operating in line with the memoryless property, if current solutions do not fare well against their counterparts in the original group, their success probabilities in subsequent iterations match those of the present.Therefore, past failures do not dictate future outcomes.Within the memoryless matrix,  www.nature.com/scientificreports/there are two categories of solutions: winners and losers.While losers can still contribute to the optimization process alongside winners, a solution is deemed victorious if its efficiency surpasses that of its counterpart in the X winners group.If a solution is updated in both the X winners and memoryless matrices, it a winner, otherwise, it a loser.Our designed exploitation model, focused on updating solutions adhering to the exponential distribution, relies on both winning and losing solutions.The updated solution, V timx+1 i is modelled around a specific involving random and adaptive parameters.This strives to locate the global peak near an efficient solution, creating a new solution for future populations.The delves deep, examining how winners and losers navigate within the search space, leveraging valuable data from winners: Additionally, the exponential rate, in relation to the mean, can be deduced: The exponential mean arises from averaging the guiding solution and the associated memoryless vector, which can be either a winner or loser.

EDO exploration potential
This section delves into the exploration aspect of the introduced algorithm.The exploration segment pinpoints potential areas within the search domain that likely house the ideal, global solution.The EDO exploratory blueprint derives from two prominent solutions, or winners, from the primary population, which adhere to the exponential distribution pattern: Updating the new solution involves with M time representing the average of all solutions sourced from the primary group.This average is computed by adding all exponential random variables from the same dimension and then dividing by the total population, denoted as N .The term c signifies a refined parameter, indicating the proportion of information drawn from the Z 1 and Z 2 vectors towards the contemporary solution and is formulated as: Here, d stands as a flexible parameter.Initially set to zero, it undergoes gradual decrement as time progresses.Here, 'time' alludes to the present moment, while 'Max_time' signifies the comprehensive duration or iterations.Both Z 1 and Z 2 are viewed as potential vectors, formed by: Furthermore, D 1 and D 2 delineate the distance between the average solution and the 'winners' randomly picked from the initial population.At the inception of the optimization journey, a notable disparity exists between the average solution and the standout performers.Yet, as the process nears its end, this gap between the prominent solutions and their corresponding variances narrows.To establish the Z 1 and Z 2 vectors, exploration is conducted around the average solution, aided by a pair of randomly chosen outstanding solutions.

Optimizing with exponential distribution optimizer (EDO)
The EDO method we are introducing follows a series of steps to thoroughly navigate the search space, aiming for the global optimum.Initially, we create a collection of solutions, randomly generated and marked by a wide range of values.The search process is depicted using various exponential distributions and as such, the location of each solution can be seen as random variables adhering to this distribution.We design a matrix without memory to mimic the absence of memory and initially, it mirrors the original group of solutions.Leveraging the exploratory and refining stages of our method, every solution starts moving closer to the global optimum over time.During the refining stage, the matrix without memory is used to store the outcomes from the prior step, irrespective of their past, allowing them to play a pivotal role in shaping the new solutions.This leads to the categorization of solutions into two groups: the successful ones and the unsuccessful ones.Additionally, we incorporate various properties of the exponential distribution, like its mean, rate and variance.The successful solution is guided by a leading solution, while the unsuccessful one follows the successful one, aiming to discover the global optimum (11) nearby.In the exploration stage, the fresh solution is influenced by two randomly chosen successful solutions from the initial group and the average solution.At the start, both the average solution and its variance are distant from the global optimum.However, the gap between the average solution and the global optimum narrows down until it reaches its lowest point during the optimization.A toggle parameter determines whether to embark on the exploration or refining stage, based on a probability where (a < 0.5).
After crafting the new solutions, each solution boundaries are verified and then they are stored in the matrix without memory.A selection strategy is employed to incorporate the new solutions from both stages into the initial group.If a new solution proves beneficial, it integrated into the primary group.By the optimization conclusion, all solutions cluster around the global optimum.In the best solution, both the mean and variance are anticipated to be minimal, while the scale parameter λ is expected to be significant.The pseudo code of single objective EDO shown in Algorithm 1.

Preliminaries of multi-objective optimization
In multi-objective optimization tasks (MOPs), there is a simultaneous effort to minimize or maximize at least two clashing objective functions.While a single-objective optimization effort zeroes in on one optimal solution with the prime objective function value, MOO presents a spectrum of optimal outcomes known as Pareto optimal solutions.An elaboration on the idea of domination and associated terminologies are illustrated in Fig. 2.

Multi-objective exponential distribution optimizer (MOEDO)
MOEDO algorithm starts with a random population of size N .the current generation is t, x t i and x t+1 i the i th individual at t and (t + 1) generation.u t+1 i the i th individual at the (t + 1) generation generated through the EDO algorithm and parent population P t .the fitness value of u t+1 i is f t+1 i and U t+1 is the set of u t+1 i .Then, we can calculate x t+1 i according to u t+1 i generated through the EDO algorithm and Information Feedback Mechanism (IFM) 72 Eq. (20).where x t k is the k th individual we chose from the t th generation, the fitness value of Begin from F 1 , all individuals in level 1 to l are added to S t = l i=1 F i and remaining members of R t are rejected illustrated in Fig. 3.If |S t | = N no other actions are required and the next generation is begun with P t+1 = S t directly.Otherwise, solutions in S t /F l are included in P t+1 and the remaining solutions N − l−1 i=0 |F i | are selected from F l according to the crowding distance (CD) mechanism, the way to select solutions is according to the CD of solutions in F l .The larger the crowding distance, the higher the probability of selection and check termination condition is met.If the termination condition is not satisfied, t = t + 1 than repeat and if it is satisfied, P t+1 is generated represent in Algorithm 2, it is then applied to generate a new population Q t+1 by EDO algorithm.Such a careful selection strategy is found to computational complexity of M-Objectives O N 2 M .MOEDO that incorporates proposed information feedback mechanism to effectively guide the search process, ensuring a balance between exploration and exploitation.This leads to improved convergence, coverage and diversity preservation, which are crucial aspects of multi-objective optimization.MOEDO algorithm does not require to set any new parameter other than the usual EDO parameters such as the population size, termination parameter and their associated parameters.The flow chart of MOEDO algorithm can be shown in Fig. 4.
Algorithm 2. Generation of t of MOEDO Algorithm with IFM procedure (20)

Results and discussion
In this section, we present the results derived from our research and experiments, which were conducted to evaluate the efficacy of the proposed method and showcase its capabilities.To ensure robust and statistically significant results, each experiment was replicated 30 times independently with population size = 40 and maximum number of iterations = 500.All techniques were executed using MATLAB R2020a on an Intel Core i7 computer equipped with a 1.80 GHz processor and 8 GB RAM.

Description of benchmark test functions
Our experimental evaluation of the MOEDO algorithm's performance utilized three separate benchmark test set functions: • Zitzler-Deb-Thiele (ZDT) test suite: from this benchmark collection 31 , we selected ZDT1, ZDT2, ZDT3, ZDT4, ZDT 5 and ZDT6 (Appendix A).A brief overview of these problems is as follows: ZDT1 presents a    www.nature.com/scientificreports/

Performance metrics for evaluation
For this research, we employed six performance metrics: hypervolume (HV), generational distance (GD), inverted generational distance (IGD), spread (SD), spacing (SP) and run time (RT).Subsequently, we provide a concise explanation of each metric to enhance comprehension shown in Fig. 5. a statistical evaluation."+/−/~" Wilcoxon signed-rank test (WSRT) was conducted at a significance level of 0.05 between the total amount of test problems on which the corresponding optimizers has a better performance, a worse performance and an equal performance for solving MO problems.

Performance analysis on benchmark test functions
Analysis using the GD metric Table 1 showcases the final solution distributions when juxtaposing our MOEDO algorithm against MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO using the GD metric.Within the ZDT benchmark, our algorithm emerged as the top performer, particularly excelling in ZDT1, ZDT2, ZDT3 and ZDT6 in terms of mean and standard deviation.In the DTLZ benchmark, MOEDO was showed commendable performance in 9/14 problems.
For the constraint test suite, the third benchmark, MOEDO outshined in both the best and average metrics, especially in SRN, OSY and KITA.A WRST test was conducted to determine the overall standing of each multiobjective method concerning the GD metric across all benchmarks.The results, presented in    www.nature.com/scientificreports/Analysis using the spacing (SP) metric Table 3 displays the final solution distributions when comparing MOEDO algorithm against MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO using the spacing metric.MOEDO showcased commendable performance in ZDT1 for average metrics, while in ZDT4, it excelled in the average metric.For DTLZ 9/14, our method stood out in the standard metric.In the constraint benchmark, encompassing TNK, OSY, BNH and KITA, MOEDO consistently scored higher across all metrics.MOEA/D lagged behind in all benchmarks.The WRST test for the Spacing metric positioned MOEDO as the top performer, followed by MOGNDO and NSGA-II.

Analysis using the spread (SD) metric
Table 4 presents the final solution distributions of MOEDO algorithm against MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO using the Spread metric.MOEDO was dominant in ZDT2 across all metrics and in ZDT2, ZDT3, ZDT4 and ZDT5 for both the std and average metrics.In the DTLZ benchmark, our method was superior in 10/14 for both std and average metrics.For the constraint benchmark, MOEDO outshined in TNK and OSY for the standard metric.MOEA/D consistently underperformed.The WRST test, based on the results from Table 4, highlighted MOEDO as the leading method among the evaluated algorithms.www.nature.com/scientificreports/Analysis using the HV metric Table 5 delineates the final solution distributions when contrasting MOEDO algorithm against MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO using the HV metric.Our method outshone in 15 out of the 26 problems.It was particularly dominant in ZDT3 and ZDT5 across metrics.In the DTLZ benchmark, MOEDO was unparalleled, especially in DTLZ4.Lastly, in the constraint test suites, MOEDO excelled in TNK, OSY and KITA across all metrics.In the ZDT benchmarks, MOAOA consistently underperformed.The WRST test outcomes for the HV metric positioned MOEDO in the second spot among the five evaluated multi-objective methods.

Analysis using the RT metric
Vol.:(0123456789) www.nature.com/scientificreports/leading multi-objective methods using metrics like GD, IGD, HV, Spacing, Spread and RT across standard test suites such as ZDT, Constraint and DTLZ.MOEDO commendable performance is attributed to its rapid convergence and a harmonious balance between exploration and exploitation.This is achieved by integrating the strengths of the IFM technique, the adept search capabilities of the EDO algorithm and the enhanced convergence facilitated by a IFM mechanism.Conversely, many of the compared multi-objective methods struggled with maintaining a balanced exploration-exploitation ratio and lacked efficient convergence.

Application to real-world engineering design problems
Table 7 presents the performance of MOEDO concerning the spacing metric, compared against MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO.A lower spacing value denotes superior coverage performance.The data reveals MOEDO performance RWMOP1, RWMOP2, RWMOP3 and RWMOP4 was suboptimal.The WRST test, as illustrated in Table 7, positions MOEDO as the top performer in terms of spacing values among all others.Table 8 showcases the findings of MOEDO for the HV metric, juxtaposed against other algorithms.A higher HV value is indicative of better performance.The data suggests that MOEDO excels across RWMOP1, RWMOP2, RWMOP5 and RWMOP6.However, the WRST test results in Table 8 indicate room for improvement in the HV metric for MOEDO.www.nature.com/scientificreports/Table 9 showcases the performance of MOEDO concerning the run time metric, juxtaposed against other MOEDO others MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO.Ideally, a lower RT metric value is preferred.Specifically, it showcased superior results for RWMOP2, RWMOP3, RWMOP4 and RWMOP5 in both Average and Best metrics.The overall ranking of each algorithm, based on the Best and Avg results, was determined using the WRST test.As illustrated in the concluding row of Table 5, MOEDO emerged as the top performer for the RT metric among the others.
Figure 10 emphasizes MOEDO superior performance for the RWMOP function, as its results closely align with the PF.

Conclusion
In this research, we introduce a multi-objective variant of the exponential distribution optimizer (EDO) method, a recently proposed metaheuristic algorithm rooted in specific principles of mathematics exponential distribution theory.We present a multi-objective EDO termed MOEDO, which integrates concepts of multi-objectivity, NDS, CD and IFM theory into the conventional EDO framework.The IFM approach, drawing from diverse strategies, ensures a balanced exploration-exploitation dynamic, promoting enhanced convergence and the ability to bypass local minima.experimental outcomes indicate that our MOEDO algorithm surpasses MOMPA,       www.nature.com/scientificreports/

Figure 2 .
Figure 2. Multi-objective all definitions in search space of MO-problem.

Figure 3 .
Figure 3.The procedure of the NDS approach based on MOEDO algorithm.

Figure 5 .
Figure 5. Mathematical and schematic view of the GD, IGD, SP, SD and HV metrics.

Figure 7 .
Figure 7. Best Pareto optimal front obtained by the MOEDO algorithm on DTLZ1-DTLZ7 problems with 2-objectives.

Figure 8 .
Figure 8. Best Pareto optimal front obtained by the MOEDO algorithm on DTLZ1-DTLZ7 problems with 3-objectives.

Figure 9 .
Figure 9. Best Pareto optimal front obtained by the MOEDO algorithm on constrained CONSTR, TANK, SRN, OSY, BIN and KITA.

Figure 10 .
Figure 10.Best Pareto optimal front obtained by the MOEDO algorithm on real-world engineering problems (RWMOP1-RWMOP5).
Table 1's, place MOEDO at the top and MOMPA bottom behind other algorithms across all benchmarks.Analysis using the IGD metricTable2presents the final solution distributions when comparing MOEDO algorithm against MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO using the IGD metric.For the ZDT3 and ZDT 6 test function, our algorithm surpassed its counterparts.In the DTLZ benchmark, MOEDO was notably superior in 8/14 metrics.Furthermore, in the third benchmark, MOEDO exhibited superior performance, especially in TNK, OSY and KITA both in the best and average metrics.MOMPA and MOEA/D consistently underperformed in all benchmarks.The WRST test's results, displayed in Table2's final row, highlight MOEDO's dominance, securing the top rank among the evaluated multi-objective variants.
Table 1.Results of GD metric of different multi-objective algorithms on ZDT, DTLZ and constraint benchmark problems.

Table 2 .
Results of IGD metric of different multi-objective algorithms on ZDT, DTLZ and constraint benchmark problems.

Table 3 .
Table 6illustrates the final solution distributions of MOEDO algorithm against MOMPA, NSGA-II, MOAOA, MOEA/D and MOGNDO using the initial gap metric.MOEDO surpassed its counterparts in DTLZ3 and DTLZ5 across all metrics, namely Best, Avg and Std.In the constraint functions, our method was superior in CONSTR, TNK, SRN and KITA for the RT metric.MOEA/D consistently trailed its peers.The WRST test, based on the results from Table6, positioned MOEDO at the best among the tested methods.The visual representation of the final solutions from our proposed MOEDO method.Figures 6, 7, 8 and 9 reveals MOEDO superior PF.This section provides a thorough assessment of MOEDO performance against Results of SP metric of different multi-objective algorithms on ZDT, DTLZ and constraint benchmark problems.

Table 4 .
Results of SD metric of different multi-objective algorithms on ZDT, DTLZ and constraint benchmark problems.

Table 5 .
Results of HV metric of different multi-objective algorithms on ZDT, DTLZ and constraint benchmark problems.

Table 6 .
Results of RT metric of different multi-objective algorithms on ZDT, DTLZ and constraint benchmark problems.

Table 8 .
Results of HV metric of different multi-objective algorithms on engineering design problems.

Table 9 .
Results of RT metric of different multi-objective algorithms on engineering design problems.